Structural Symmetry of Two-dimensional Metallic Arrays: Implications for Surface 

Plasmon Excitations 
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In recent years, there has been intensive investigation of surface plasmon polaritons (SPPs) in the 
science and engineering fields. Understanding the physics of surface plasmon excitation is essential 
to the manipulation of SPPs, and most existing studies focus on (— l,0)-type SPP excitation. In 
this article, we report our recent investigation of the (0, ±l)-type SPP excitation of a gold two- 
dimensional nano-cavity array using finite- difference time-domain methodology. Our particular 
focus is on the symmetry properties of (0, ±l)-SPPs excited by different polarizations of light. It is 
found that polarization has strong implications for the field distribution of the corresponding SPPs. 
As a result, the control of polarization may provide important insights into the manipulation of 
SPPs. 
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1. INTRODUCTION 

Recent advances in nanotechnology have ignited the 
study of surface plasmon polaritons (SPPs), which have 
become a popular topic [3, 2, 3] due to their wide range of 
potential applications in surface-enhanced Raman scat- 
tering (SERS) [J], surface enhanced second harmonic 
generation, nano-photonics [5[, thermal-photovoltaic de- 
vices 0, 0, H, Q , and data storage and imaging [Tol ] 

However, finding a way to excite and control SPPs in 
an advantageous manner is still the main concern of sci- 
entists. The technology available to control SPPs with 
precision and flexibility remains underdeveloped, and dif- 
ferent schemes have been proposed in this regard. For ex- 
ample, since the discovery of extraordinary transmission 
major interest has centered on investigating cylindri- 
cal hole arrays in which subwavelength hollow cylinders 
are periodically formed on a fiat metal film using litho- 
graphic methods 0, El- The work of Keif et al. pj|, 
however, has demonstrated showed that the shape of an 
individual cavity in metallic grating plays a dominant 
role in controlling the excitation of SPPs. In addition, 
van der Molen et al. [15] showed that there are shape 
and localized resonances in two-dimensional (2D) peri- 
odic subwavelength metallic cavity arrays. As a result, 
different shapes can lead to different resonances because 
the holes act like plasmonic cavities to confine the electric 
field and thereby give rise to strong localized resonances 
01 EH E3 • Understanding the geometry effects of cavi- 
ties can lead to a number of applications that require the 
precise spatial and frequency control of an enhanced elec- 
tric field, such as SERS and thermal radiation. However, 
scientists have only a limited understanding of the sur- 
face shape resonances, and the correlation between SPP 
excitation and cavity geometry has not been widely in- 



vestigated until recently [3, El, El E3, IH S [H, \jZ^. 
In a recent work, we reported the fabrication of 2D nano- 
cavity arrays on a gold surface using interference lithog- 
raphy (IL) [2l|. In this article, we focus on a theoretical 
study of 2D nano-cavity arrays by using finite- difference 
time-domain (FDTD) simulation [25[ methodology. Our 
particular focus is on the symmetric properties and po- 
larization dependence of excited SPPs 



2. SIMULATION DETAILS 



2.1. Basic simulation cell setup 
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FIG. 1: (Color online) 
cylindrical cavity. 



a) Unit cell and (b) top view of the 
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We consider a model system that contains a gold slab 
drilled with a 2D array of cylindrical cavities. It is con- 
venient to define the unit cell shown in Fig. [1] with 
suitable boundary conditions. Periodic boundary con- 
ditions (PBCs) are applied in both the d=x- and in- 
directions to produce a periodic structure. To account for 
non-periodici ty in the z-direction, the perfectly matched 
layer (PML) [26[ boundary condition is adopted in the 
indirections. In addition, certain parameters are fixed 
throughout the study. For example, we only examine 
optical wavelengths that range from 400 nm to 900 nm. 
Moreover, the thickness of the gold layer is set at 1 jiva 
to ensure that the structure can be considered to be op- 
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tically thick and such bulk photonic effects as Fabry- 
Perot resonance and guiding mode resonance [27[ are 
substantially reduced and negligible. Finally, the spatial 
and time resolution are Ax = 25 nm and At = Ax/c, 
respectively, and the two lattice constants are fixed at 
a = 575 nm. 
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2.2. Dielectric function of gold 



TABLE I: (Color online) Fitted parameters for Eq. 0] 



Gold processes complex electron inter-band transitions 
at infrared frequencies. As a result, simple dielectric 
functions such as the Drude model, which considers only 
the electric response of conduction electrons, may not be 
adequate for modeling the dielectric response of gold. In 
this work, the multiple Drude-Lorentz model proposed 
by Vial et al. [28[ is adopted: 



(1) 



Here, index j represents the contribution to the dielec- 
tric function made by the electrons of different bands; ujj 
and 7j are the frequency and damping parameters; 
is the high-frequency response, which is originated by 
the screening of the core electrons; and, finally, Asj = 
e oo — £jj where Sj is the static dielectric response con- 
tributed by the different electron bands. Equation (pQ) 
is fitted against the measured dielectric function of gold 
[29[ , and the fitted parameters are tabulated in Table H 
Figure [2] shows both the measured and fitted dielectric 
functions. As can been seen, Eq. (pQ) agrees well with the 
measured e throughout our frequency range of interest. 
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FIG. 2: (Color online) The real (e r ) and imaginary (si) parts 
of the dielectric function fitted with the multiple Drude- 
Lorentz model. 



entire optical frequency range. They are randomly lo- 
cated so that all of the resonant modes of the system 
can be excited. After turning off the sources, some fields 
are left in the system to allow their magnitudes to be 
recorded as a finite-length signal. This signal is then ex- 
pressed as the sum of a finite number of sinusoids with 
exponentially decaying factors in a given bandwidth. The 
frequencies, decay constants, amplitudes, and phases of 
these sinusoids are determined from the coefficients of the 
sum. 

Furthermore, line sources pointing in the ^-direction 
are used to allow an investigation of the effects of dif- 
ferent polarizations on the dispersion relations. These 
sources can be classified as p-polarized or s-polarized, 
depending on whether it is the magnetic field (H y ) or 
electric field (E y ) that pointis in the ^/-direction, respec- 
tively. Some of the resonant modes shown in Fig.[3]corre- 





0.05 0.1 0.15 0.2 0.25 
k x (27t/«) 



FIG. 3: (Color online) Dispersion relations of an r = 75 
nm cylindrical cavity array under (a) p-polarized and (b) s- 
polarized lights. The green dotted lines are Wood's anomalies 
[Eq. J2}], the red lines are the SPP dispersions of a gold slab 
[Eq. ©], and the black dotted lines are the calculated (0, ±1) 
modes. 

spond to Wood's anomalies [l2[. Wood's anomalies occur 
when one diffracted order is parallel to the structure sur- 
face, which results in energy redistribution to the other 
diffracted orders. In other words, the incident light is 
scattered in parallel to the structure surface. This causes 
a sudden change in reflectance when the frequency (uj) 
and in-plane wave vector (k x ) of the incident light satis- 
fies 



2.3. Dispersion relation calculation 



uo = c\k\ = c\h x + nG x + rriG 



2/h 



(2) 



To calculate the dispersion relation, several Gaussian 
sources are placed within the simulation cell. These 
Gaussian sources have a frequency width that covers the 



where k is the incident wave- vector, G x and G y are the 
reciprocal lattice vectors, and n and m are integers. Ad- 
ditionally, we can trace the SPPs by plotting their dis- 
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persion relation driven by the following equation 



UO = c 



£Au 



^Au^air 



-\k x + nG x + mG y \. 



(3) 



It can be seen that Eq. (j3]) agrees well with the remain- 
ing excited modes of the p-excited resonances for (n, m) = 
(— 1, 0), (1, 0), and (0, ±1), whereas only the (0, ±l)-type 
resonances are present for s-polarization. This is due to 
the fact that the array has C± v point group symmetry. 
The eigenstates with k parallel to the Y-X direction are 
irreducible representations of point group C\h [30], and 
each eigenstate is either symmetric (A mode) or antisym- 
metric (B mode) under a mirror transformation along the 
x-z plane and can be excited only by incident waves that 
have the same symmetric property. These are, as is dis- 
cussed in Appendix O p- and s-polarized lights, respec- 
tively. Furthermore, it is proven that only a symmetric 
eigenstate is present for (n, m) = (—1,0) and (1,0), but 
that, on the contrary, both symmetric and antisymmet- 
ric eigenstates exist for (n, m) = (0,±1). Therefore, the 
black dotted lines shown in Figs. 02(a) and (b) represent 
the eigenstates of the A and B symmetries, respectively. 



3. RESULTS AND DISCUSSION 



3.1. Field pattern calculation 



Consider the time-average of energy density, as follow 
J(E*(r,t)-E(r,t))dt 
= J (J E*(r,uj)e- iujt duj • J E(r,uo')e iuj/t duo')dt 
(E*(r,u) ■ E(r,u'))( J e^' dt)duodJ 
{E*(r,u) ■ E{r,uo f ))2ir5(uo' - uo)dwduo' 
2tt / (£*(r>) • E(r,uo))duo. 



Spectral field density is a single uo component of & (r, uo) • 
E(r,uo). Theoretically, we need to sum up all of the uo 
components. However, as the Gaussian source used has 
only a narrow frequency width, Acj, the component of the 
specified central frequency, uj cen dominates the integral. 
By assuming the other ethers components are weak, that 
is, 



\E(r,uo cen )\ 2 >> \E(r, Others) | 2 , 



(5) 



we have 



|£(r,^ cen )| 2 = J E*(r : t).E(r,t)dt. (6) 

In addition, \E(r, oj cen )\ 2 also provides information about 
the harmonic eigenstate with a specified k x and cc; cen . 



We are now in a position to present the field density 
distribution of the various excited SPP modes. Guided 
by the dispersion relations, we have located the wave- 
vectors and frequencies of these excited modes at k x = 
0.1725, uo = 0.9968 for the p-excited resonance from Fig. 
[3](a) and k x = 0.1725, uo = 0.9849 for the s-excited res- 
onance from Fig. [3^b) , respectively, and have calcu- 
lated their corresponding near field intensities. To cal- 
culate field intensity, a plane-wave Gaussian source is 
put above the array at z = 1 /im to mimic an incident 
beam. The source has a narrow frequency width, such 
that it simulates a single- frequency incident plane wave. 
In addition, a spatially dependent amplitude function, 
ex.-p(27rik x x), is added to the source, where k x is the in- 
plane wave- vector, such that the angle of incidence is 
given by sin -1 (k x c/uo). 

This simulation provides us with useful information 
on the way in which the incident light interacts with 
the array and the physical properties of the excited res- 
onances. Time-domain simulations provide information 
about both the transient and steady states. We focus on 
the steady states by calculating the spectral field density 
of an excited mode at specified k Xl uo values: 

\E(f,oj)\ 2 = [ (E*(f,t) ■ E(r,t))e^dt. (4) 



3.2. Symmetry properties of the excited SPPs 

In this section, we present the calculated field densities 
of the excited (0, ±l)-SPPs and discuss their symmetry 
properties. A general discussion of the symmetric prop- 
erties of the eigenstates [30|, [HI is given in Appendix S 
Electric field density (\E(r, uo)\ 2 ) under p- and s-polarized 
excitations is shown in Figs. Hfa)-(d) and Figs. E^a)-(d), 
respectively. The colors indicate field strength, which 
ranges from 0-80 arbitrary units (A.U.). The cavity is 
outlined with black dotted lines. As can been seen in 
Fig. Ufa), strong electric fields are distributed along the 
x-axis (y = 0), with the strongest being concentrated 
near the rim of the cavity. This is in accordance with 
Eqs. (lMB - (lAT3l) . where E 2 X [Fig. [1(b)] and E 2 Z [Fig. 
21(d)] are symmetric along the x-z plane, and E 2 [Fig. 
21(c)] is anti- symmetric. Therefore, E 2 — along y = 
and y = ±a/2. In contrast, the field is asymmetric along 
the x = plane, which is due to the propagating na- 
ture of the excited SPP mode [22[ . When calculating the 
\E(r, oo)\ 2 , two kinds of electric fields are present in the 
system: one is the electric field of the resonant mode, 
labeled E res (r,uo), and the other is the electric field of 
the incident light, labeled E[ nc (r, uo). According to group 
theory, these two electric fields should have the same 
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FIG. 4: (Color online) Electric field density on z — of the 
(0, ±1) SPP mode of an r = 75 nm cavity array excited by a p- 
polarized light: (a) \E(r,u)\ 2 , (b) \E x (f,cv)\ 2 , (c) \E y (r,u)\ 2 , 
and (d) \E z (f, oj) | 2 . The in-plane wavevector and frequency 
of the excited SPP are 0.1275 and 0.9968, respectively. The 
fields are not drawn to scale for better visibility. Note the 
symmetry along y = and y — zba/2. 



be dominated by point dipole-like distribution, (that is, 
E x oc p x ), where p x is the total dipole moment along 
the x-axis. This is in accordance with the hypothesis of 
Teperik et al. [32], who treats each cavity as a resonator 
that consists of a simple capacitor, an inductor and re- 
sistance. 
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FIG. 5: (Color online) Snapshots of the electric field distri- 
butions, E(r, £), for the p-excited SPP mode of Fig. |4] shown 
in (a) and (b), which are E x (r,t) and E y (r,t), respectively 
(E y is weak and can barely be observed). Note the symmetry 
along y = 0. Blue and red represent positive and negative 
values, respectively. 



FIG. 6: (Color online) Electric field density on z — of the 
s-excited (0 ± 1) SPP mode, k x = 0.1275 and u = 0.9849. 
The symmetry along x = and y = indicates the localized 
nature of the resonance. 





(b) 



symmetric properties under C\h symmetric operations. 
The p-polarized light contains an electric field at the in- 
cident plane (the x-z plane) whose time average is always 
symmetric along the y = plane, but asymmetric along 
the x = plane. As a result, the resultant field density, 
\E(r, oj) | 2 , is asymmetric along the x = plane. 

To further investigate the p-excited (0, ±1) SPP, we 
show in Fig. [5] a snapshot of the field: E x (r,i) and 
E y (r,t) during the FDTD simulations. The symmetry 
conditions, E x (x,y,z) = E x (x,—y,z) [Eq. (|Alip ] and 
E y {x 1 y 1 z) = —E y (x,—y,z) [Eq. (|A12|) ]. can be ob- 
served clearly in the figures, although E y is too week 
and can barely be observed. It should be noted that the 
x-component of the electric field is very weak at x = 0, 
and is mainly concentrated at (x = ±r), which suggests 
that the local field should, at least to the first order, 



FIG. 7: (Color online) Snapshots of the electric field distri- 
bution of the s-excited SPP mode in Fig. [6] 



We next turn our attention to the s-excited (0, ±1) 
SPP modes, \E(r, cj)| 2 , under s-polarized excitation, as 
shown in Figs. IEKa)-(d). According to Eqs. ()A11|) — 
(|A13|) . E x and E z should be anti-symmetric and E y 
should be symmetric along y = 0, which can be ob- 
served in Figs. E3(b)-(d). It should be noted that the 
field densities of E% and Ey are very similar. However, 
as is revealed in Figs. H(a)-(b), E x and E y are actually 
very different from each other. Moreover, it should also 
be noted that E y can be attributed to an effective dipole 
along the ?/-axis. In other words, E y (xp yi which is simi- 
lar to the case for E x in the p-excited (0, ±1) mode [Figs. 
SKa) and[5ja)]. 
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3.3. Conclusion 

The present analysis relies only on the knowledge of 
C^-point group symmetry, which should be valid for 
any plasmonic structure that belongs to the same point 
group. Therefore, our argument can be applied to all 
C^v -point group structures, such as square cavities, cross- 
shape grooves, circular crafters, pyramidal structures 
[20l . l24j , and bottle-shape cavities [2]], [22[ . Additionally, 
SPP excitation along the T-M direction also has similar 
symmetric properties [30], and the corresponding eigen- 
modes can also be classified into even and odd modes, 
which are excited by p- and s-polarized light accordingly. 

To sum up, we have performed FDTD analysis on 
a gold cylindrical cavity array and investigated the ex- 
cited (0, ±1) mode based on the symmetry argument. 
Although very close in resonance frequencies, the even 
and odd modes display very different types of behavior 
in terms of field distribution, which has important im- 
plications for near-field applications such as SERS. It is 
hoped that this work will provide insight into the fields 
of plasmonics and nano-optics. 



APPENDIX A: SYMMETRIC OPERATIONS OF 
EIGENMODES 

The eigenfunctions, H^(f,u), of a periodic plasmonic 
structure at eigenfrequency uj are Bloch functions that 
satisfy the master equation [50, [3l| : 

ei?£(r» = Vx(lvx)%» 



However, vectors defined by cross products are trans- 
formed differently. For example, for magnetic field 



(Al) 



where k denotes wavevectors within the first Brillouin 
zone. Consider a mirror transformation operator, M xz , 
along the x-z plane that is mapped from R 3 to R 3 by 

X Z X X , 

V M xz y -y, 
z — M xz z — z. 

Under this transformation, the ^/-coordinate obtains a 
minus sign, and the x- and ^-coordinates remain the 
same. The vector fields are transformed by operator O xz , 
defined as 

O xz (f(r)F(r)) = f(M xz r)M xz (F(M xz r)), (A2) 
and thus, 

O xz f(r) = f(x,-y,z), (A3) 
O xz E x (x,y,z) = E x (x,-y,z), (A4) 
O xz E y (x,y,z) = -E y (x,-y,z), (A5) 
O xz E z (x,y,z) = E z (x,-y,z). (A6) 



F x x + F y y 



F = qv x B, 
F z z = 



x y z 

V X Vy V Z 

B X By B z 



(v y B z - v z B y )x 
+ (v z B x - v x B z )y 

+ (V X By ~ VyB X )Z. 

To preserve the transformation properties of the force 
and velocity, the magnetic field is transformed as follows 



O xz B x (x,y,z) 
O xz B y (x,y,z) 
O xz B z (x,y,z) 



-B x {x,-y,z), (A7) 
B y (x,-y,z), (A8) 
-B z (x,-y,z). (A9) 



Electric fields and magnetic fields are transformed in 
different ways under the mirror transformation. 

It is well known that commutes with O xz for periodic 
structures: 



eo xz = o xz e. 

Also, if k — kx, then we have 

O xz k = /c, 

and 

6 xz H^{x, y, z) = ±H%(x, y, z). 



(A10) 

That is, the eigenmode is either symmetric (+) or an- 
tisymmetric (— ) under mirror transformation along the 
x-z plane. The substitution of Eqs. (|A4|) - (|A6|) and (|A7|) - 
(lA9l) into Eq. (|AlO|) gives 



El 


:X ( x ,y, z ) 




-y,z), 


(AH) 


E l 


, y (x,y,z) 


= 


-y,z), 


(A12) 


E l 


>z ( x ,y, z ) 




-y,z), 


(A13) 


H % 




= TH s Jx, 


-y,z), 


(A14) 


H% 


, y (x,y,z) 


= ±Hr (X, 

k,y v ' 


-y,z), 


(A15) 


Hi 


z ( x ,y,z) 


= TH % Jx, 


-y,z) 


(A16) 



for eigenmodes that satisfy Eq. (|A10|) . In addition, for 
p-polarized lights, E = (E x ,0,E z ), H = (0,i^,0), and 

E x (x,y,z) = E x (x,-y,z), 
E z (x,y,z) = E z (x,-y,z), 
H y (x,y,z) = H y (x,-y,z). (A17) 

In contrast, E = (0,£^,0) and H = (H x ,0, H z ) for s- 
polarized lights; therefore 

E y (x,y,z) = E y (x,-y,z), 
H x (x,y,z) = H x (x,-y,z), 
H z (x,y,z) = H z (x,-y,z). (A18) 
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A com parison of Eqs. (|A17j) and (|A18j) to Eqs. (|Allj) - 
(|A16j) suggests the p- and s-polarized lights are sym- 
metric and antisymmetric to the x-z plane mirror trans- 
formation, respectively. Therefore, the symmetric (+) 
modes can be excited only by p-polarized light that con- 
tains a symmetric electric field along the x-z plane. Sim- 
ilarly, antisymmetric modes (— ) can be excited only by 
5-polarized light. 
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